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abstract 


Pre-launch and post-launch calibration plans for the APAFO magnetometer 
experiment are presented. A study of trade-offs between boom length and 
spacecraft field is described; the results are summarized. The , 5 ™^ star 
plan includes: calibration of the Vector Fluxgate Magnetometer (VFM) Star 
Sensors, and Scalar Helium Magnetometer (SHM); optical bench integrati , 
and acquisition of basic spacecraft field data Post-launch calibration 
has two phases. In phase one, SHM data are used to calibrate the VFM total 
vector magnetic-field data are used to calibrate a physical model of the 
spacecraft field, and both calibrations are refined by iteration. In phase 
two corrected vector data are transformed into geocentric coordinates, 
previously undetected spacecraft fields are isolated, and initial geomagnetic 
field model, are computed. Provided the SHM i, accurate to the required 1.0 nT 
and can be used to calibrate the VFM to the required 3.0 nT accuracy the 
trade-off study indicates that a 12-m boom and a spacecraft field model 
uncertainty of 5Z together allow the 1.0-nT spacecraft field error 
requirement to be met. 
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INTRODUCTION 


The APAFO magnetometer experiment consists of a three-axis vector 
fluxgate magnetometer and a scalar helium magnetometer together with two 
star sensors. These are mounted on an optical bench at the end of a boom 
in order to reduce any magnetic fields from the spacecraft, its systems and 
other investigations. At present, the length of the boom is baselined at 
12 m (see Figure 1). The scalar magnetometer sensor will be provided by E. 
Smith of the Jet Propulsion Laboratory (JPL); the vector magnetometer 
sensor and the electronics for both magnetometers will be provided by M. 
Acuna of Goddard Space Flight Center (GSFC); the star sensors will be 
provided by T. Potemra of The Johns Hopkins University Applied Physics 
Laboratory (APL) . All three are APAFO Co-Investigators. The optical bench 
and boom will be provided by vendors to be selected during Phase C. 

The error budget for the magnetometer experiment, in nanoTeslas (nT) , is 
as follows: 


SCALAR VECTOR 



Req. 

Goal 

Req. 

Goal 

Instrument 

1.0 

0.7 

3.0 

1.5 

Spacecraft Field 

1.0 

0.5 

1.0 

0.5 

Position and Time 

1.3 

0.5 

2.0 

1.0 

APAFO Field Error 

0.2 

0.1 

0.2 

0.1 

Attitude Error 

— 

” ~ ~ 

3.0 

2.0 

Rss 

1793 

T7o 

O 

1774 


It is our intention to meet the goal and not merely the requirement. To do 
so in practice we must attempt to make the errors less than the goal. 

Toward that end, we suppose for now that the APAFO field error, i.e. that 
error due to stray magnetic fields from the instruments mounted at the 
optical bench, can be reduced to 0.1 nT or less. We may then proceed to 
consider the other sources of error. 

The purview of this report, then, are the instrument and spacecraft field 

contributions to the error budget. For this mission, the basic in-flight 
standard of accuracy is the scalar helium magnetometer (SHM) . It is 
believed that the SHM can be built and calibrated pre-launch such that its 
basic accuracy is within 0.7 nT . This report shows that if this is 
accomplished then the meeting of the error budget goal for spacecraft field 
and APAFO field is feasible, provided an adequate model of the spacecraft 
field is derived and provided that the input information needed in that 
model is present in the spacecraft telemetry. 

To accomplish this will require certain steps be taken both prior to 
launch and in orbit to calibrate both the magnetometers and the spacecraft 
field model. These procedures are the topic of this report. 
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PRE-LAUNCH CALIBRATION 


Vector Fluxgate Magnetometer (VFM) 

1. Internal calibration parameters will be established by comparison 
with an absolute instrument at the GSFC magnetic test site. 

2. The relationship of the internal magnetometer axes to a coordinate 
system determined by an externally mounted optical cube must be 
experimentally determined. 

APAFO Magnetic Fields 

An attempt will be made to determine parameters for fields exceeding 0.1 
nT at the SHM due to the VFM and star sensors and for fields exceeding 0.1 
nT at the VFM due to the star sensors. 

Star Sensors 


APL and their external vendors shall furnish calibration parameters 
relating the field of view of the sensor to a coordinate system determined 
by an externally mounted optical cube. 


Scalar Helium Magnetometer (SHM) 

The accuracy of the SHM shall be determined by laboratory measurement of 
the various parameters contributing to instrument error. This error budget 
is made up of: 

Light Shift 0.5 nT 

Noise 0.15 nT 

Servo Error 0.1 nT 

Readout Error 0.07 nT 

Electronic Offset 0.05 nT 

which gives a maximum error of 0.87 nT or an rss error of 0.54 nT. 

In addition, the overall error shall be verified by direct comparison 
with a Proton Precession Magnetometer over its range of operation and over 
a comprehensive selection of field directions with respect to the 
magnetometer axis. 


Integrated Optical Bench 

Optical alignments shall be determined between the four optical cubes: 
one each on the two star sensors, one on the VFM and one on the SHM. This 
should be accomplished in the GSFC optical facility according to a plan 
submitted by the contractor responsible for integration and test. 

Optical alignment accuracies are expected to be within 1 arc-second, 3a. 
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Spacecraft 


1. The locations of the magnetic torquers, solar panels, and any other 
magnetic sources likely to cause a field >0.1 nT at the optical 
bench, shall be measured relative to the base of the APAFO boom to an 
accuracy of about 3-5 cm. 

2. The pole-strength of the magnetic torquers as a function of the 
current applied shall be measured to at least 1 Z. 

3. As possible, the magnetic characteristics of other magnetic sources, 
including investigator instruments, shall be specified (measured) as 
a function of their operational state. 

4. The housekeeping telemetry shall include, as possible, 

a) Torquer current, 

b) Solar panel current, position and other state 
information, 

c) Main bus current, 

d) Pertinent state information from other identified 
magnetic sources. 

5. The boom will need to be deployed to a position whose accuracy should 
be known to about 2-3 cm. 
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POST-LAUNCH CALIBRATION 


It cannot be assumed that the pre-launch APAFO magnetometer calibration 
will be unchanged by launch, boom deployment, and long-duration exposure to 
the space environment. In-flight magnetometer calibration is necessary to 
monitor and correct for changes during launch and on orbit. Post- launch 
calibration will take place periodically during the mission. The frequency 
of the various steps in this calibration will be determined after launch 
when it is found how rapidly the calibration parameters are changing. 
Calibration must enable accurate calculation of the geophysical parameters 
(the ambient magnetic flux density vector) from APAFO data both during and 
between calibrations despite variable spacecraft fields, small boom 
motions, and VFM drift. Periodic in-flight calibration of the VFM against 
the SHM is used to derive the VFM calibration parameters. These 
periodically updated parameters enable calculation of total vector magnetic- 
field values between calibration runs; however, calculation of the ambient 
field also requires elimination of the spacecraft field from the observed 
total field. This requires a model of the spacecraft field which itself 
must be calibrated. 

Post-launch calibrations will consist of two phases, each with several 
steps. The first phase consists of those steps necessary to determine the 
VFM calibration parameters and the spacecraft field model parameters. The 
second phase consists of those steps necessary to confirm and/or modify 
transformations between spacecraft and Earth-centered coordinates and to 
model any previously unaccounted-for fields originating on the spacecraft. 


PHASE ONE 


At present this phase will be described as a multi-step iterative 
process. It is easier to comprehend and describe when viewed in this 
manner. However, in practice the procedure may be integrated. The actual 
procedure to be followed will be determined in a Phase C study. 

Phase One: Notation and Coordinates 


Begin by introducing time, t, dependent coordinate systems and position 
vectors therein: 

r(t) = position of the ESA/PPF spacecraft coordinate system origin in 
Earth-centered spherical polar coordinates. 

x°(t) = position of the optical bench coordinate system origin in 
the spacecraft coordinate system (x--12m,y-8m, z-Om) . 

xj (t) = position of the j th sensor in spacecraft coordinates. 

x*J(t) = position of the sensor in optical bench coordinates. 

x*(t) = nominal position of the observation in optical bench 

coordinates, the ’center’ of the SHM (x* = x*^). 
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Assume a rigid optical bench: x*J(t) - x*3 and x*(t) - x . Boom flex, 
twist, and stretch indicated by time-varying x°(t) admit time-varying xJ(t) 
despite steady x*J . 

The true magnetic flux density vector B is the sum of the ambient field 
rA which would exist in the absence of the spacecraft, the spacecraft field 
B s , and the platform or bench field B p generated by APAFO instrumentation 
and equipment: 

A S P _ 

B(x*J,t) = B (x*J , t) + B (x*J,t) + B (x*J,t) . (0) 

The transformation of (0) between coordinate systems is a tri-axial, 
time-dependent rotation. The default coordinate system is the optical 
bench coordinate system, which must be transformed into a celestial 
reference frame using bench-mounted star cameras and then into a 
terrestrial reference frame using ESA/PPF position data. For APAFO we 
consider, for generality, one quadraxial VFM (j = 1,2,3, 4) and a single 
cell SHM magnetometer (j = 5) on the bench. Note that quadraxial rather than 
triaxial fluxgates would be appropriate to the 4 -coil/ 2 -ring core, null- 
in-plane design. Each VFM coil is considered as defining one sensor axis. 
Two sensor axes (j = 1 and j — 4) will be very nearly parallel. Onboard or 
inboard magnetometers (j > 5) could improve spacecraft field estimation and 
elimination, but are not expected. 

Assume digital data is received from functioning magnetometers in the 
form of bit streams: 

ujk(tf) 5 kth bit of the offset field generator code for the j th 
VFM sensor at time ti (j S 4; k i 7); these bits define 
the integer gj(ti):(-64 i gj (ti) i 63). 

f i (tf) = fine measurement from the j tb VFM sensor fine reading 
- a 16- bit integer -32768 i fj i 32767 (j ^ 5) . 

f 5(ti ) = fine reading from the SHM (a 19- or 21-bit positive integer). 

I r (tf) = signed integer representing measured current 
energizing the r^h torque rod (r = 1 » 2 ). 

I 4 (ti) = positive integer representing current from the solar array. 

C + aC = He scalar calibration coefficient + uncertainty therein. 

afj(ti) = precision of fj(ti) (+ 0.5 in the least significant bit if fj 
represents a rounding to the nearest integer; truncation would 
merely contribute to VFM biases). 

Ol r = precision of Ir? an d 0^4 = precision of I 4 . 

Note the O prefix indicates the uncertainty in the ensuing quantity. 

The goal is to use these data, which by design must represent the 
simultaneously measured field at times t^ (i = 1,2,3,...I) during the 
calibration run t 0 i ti i tf, to calibrate the magnetometer system. 
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Phase One: Method 


The present method extends the work of Lancaster et al., "Magsat Vector 
Magnetometer Calibration Using Magsat Geomagnetic Field Measurements" Nov. 
1980; NASA Technical Memorandum 82046. The main problems considered are 
how to: 

(I) Calibrate the VFM (j=l,2 ( 3 f 4) against the SHM (j = 5). 

(II) Calibrate the spacecraft field model so as to obtain initial estimates 
of B s (x*,t) and B^(x*,t) both during and between calibration runs. 

(III) Correct and improve (I) and (II) by iterative refinement. 


At present, only the magnetic torquers are considered primary spacecraft 
field sources. In the initial iteration of Step (I), pre-launch 
calibration parameters of the VFM will serve as the starting, or a-priori, 
model. Also, pre-launch calibration parameters for the spacecraft field 
sources may be utilized in Step (I) for the calculation of field gradients 
between the SHM and the VFM. 

The present assumptions in Step (I) are: (1) that the SHM accuracy is t 
0.7 nT or better for its entire operating range; (2) that the VFM can be 
calibrated, using the SHM, to an accuracy of 1 nT, la, in the magnitude of 
its measurement and to an accuracy of 2 nT, la, in the accuracy of its 
measurement of each component. The procedure for accomplishing this is 
described in this section and is based upon a similar procedure used for 
Magsat as described by Lancaster et al. (1980). Results of applying the 
procedure to the Magsat data are described by Langel et al. (1981). 

As will be pointed out, full calibration of the VFM (Step (I)) requires 
that, during the collection of data for a calibration run, the ambient 
field along each VFM axis swing through its full range while data are 
available from the scalar instrument. Such full-range calibration is 
highly desirable but not essential. What is essential is that the 
appropriate combinations of VFM calibration parameters be determined such 
that all VFM states encountered during normal measurement of the Earth’s 
field will be calibrated. This should be achievable with appropriate 
geometric orientation of the VFM and SHM axes. Such locations will be 
determined in Phase C of the project. 

The more ambitious goal of full-range VFM calibration should also be 
pursued. Such calibration may require a planned attitude maneuver of the 
spacecraft for a short period of time (say, about 10 minutes) upon rare 
(say, once every 6 months) occasions during the mission. Attitude 
maneuvers needed to maintain orbit attitude and ground-track repeat cycle 
may prove adequate for this purpose. 
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Problem (I) is solved with estimates of steady calibration parameters, 
p n , and the time-varying total vector field values at the bench, B(x*,ti), 
with components Bik (k - 1,2,3). These estimates are to be obtained by 
iterating an optionally constrained, weighted, linearized least -squares 
calibration algorithm on data acquired at I discrete times ti during each 
calibration interval t Q i t± i tf. VFM drift is identified as any change 
in the p n which contributes more than about 0.2 nT change in measurement 
value between calibration runs. Problem (II) is solved with estimates of 
primary spacecraft field calibration parameters Mi r . These estimates may 
be derivable from routinely collected data but should be verified during 
times when a distinct, rapidly varying torque rod current signal can be 
separated from the more slowly varying ambient field. Problem (III) is 
solved by using results from (II) to correct (I). 

An advantage of this method is that VFM calibration need not depend upon 
either the spacecraft field parameterization or the twist, flex, or stretch 
of the boom. This could be achieved by solving (I) alone with data 
corresponding to times of weak or steady spacecraft field (e.g., in shadow 
with all torque rods off) . Each step also poses a smaller computational 
problem than would a simultaneous approach. The main disadvantages are as 
follows. (a) During iterated VFM calibration runs, effects of field 
gradients across the VFM caused by the spacecraft and bench instrumentation 
are treated as due in part to a spacecraft field model; the second and 
higher VFM calibration iterations thus depend weakly upon spacecraft field 
parameterization. (b) A steady field gradient appears indistinguishable 
from VFM bias. (c) Unmodeled boom motion may degrade spacecraft field 
calibration and thus, ambient field estimation. 

Problem (a) influences admissible ti, boom length, and bench fields, and 
the utility of the method. Problem (b) is addressed by the Phase two 
calibration. As to problem (c), simultaneous attitudes for both the APAFO 
bench and the spacecraft derived from star cameras mounted on the bench and 
onboard the spacecraft are to give their relative orientations to about 121 
arc-seconds rss;with a mechanical model of the boom, this relative 
orientation might be used to estimate boom-end position and remove the 
effects of translation (as well as rotation) through the spacecraft field. 
Such estimation could be difficult due to complex forcing and boom 
overtones; the boom should be designed to keep such effects small. The 
accuracy of the spacecraft field corrections, and the utility of VFM data 
to calibrate a spacecraft field model, are here considered limited by 
effects of boom-tip t ranslat ion--as distinct from the relative orientation 
of bench and spacecraft coordinate systems. Fortunately, knowledge of the 
latter allows tagging of unexpected, off-nominal, boom-tip positions. 

The method seems well suited to a long-boom, far-field case in which 
secondary spacecraft field gradients across the bench are weak and effects 
of VFM translation; i.e.,from boom motions, through the spacecraft field 
are usually small. The major difficulty is the potential inadequacy of the 
spacecraft field parameterization due to insufficient ground-based 
spacecraft field measurements, unmodeled sources, or uncertainty in the 
deployed position of the bench relative to the spacecraft. 
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Phase One: Step I, Vector FluxRate Calibration (VFC) 


This section extends the work of Lancaster et al., 1980 on the MAGSAT 
calibration. From the He scalar magnetometer 

lB(x* 5 ,ti)l = [c + ac][f 5 (ti) + 5 ( ti) 3 = c [ f 5 ( ti) + a5(ti)] (l) 

where C ^5 is the root -sum-square uncertainty in the scalar field. 

At any given time the measured field value along the jt * 1 sensor 
axis is, for j £ 4, assumed to be of the form 

K _ L 

yj = E wjkUjk + E w*j]Nji + aj(fj - bj ) (2) 

k=l 1=1 

where 

Ujk = -ujk for k £ 6 and Uj 7 = uj 7 for K = 7. 

wjk = calibration coefficient for the k th offset generator bit of the 
j th sensor coil. 

Nji = Nji(ujk) is the 1 th non-linear function representing the 

non-linear response of the jth sensor coil. The form of the 
Nji(Hjk) must be specified by ground calibration. 

w*ji = calibration coefficient for Nji — for the 1 th non-linear response 
functional of the jth sensor coil. We anticipate w* j ]Nj 1 = 0 j 1 . 

aj = scale factor for the jth sensor fine reading. 

bj = bias of the jth sensor. 

It is convenient to introduce the generic measured form 

M 

yj = a j (fj ” b j ) + E c jmf’jm( u j ) (j ^ (3a) 

m= 1 

ys = cf 5 (3b) 

where: a j , bj f j , I 5 , and C are as before; uj is the K element vector of 
bits from the offset field generator code of the jth sensor; the cj m are 
the offset field calibration coefficients (wjk or w*ji);_and Fj m (uj) is the 
mth known calibration function for the jth sensor (just Ujk if L = 0 and M 
= K) . For MAGSAT, L = 2 and Nji and Nj 2 were respectively the third and 
fifth power of signed gj (uj ) . For APAFO , we expect K = 7, guess L = 1, and 
anticipate fj to be a 16-bit result of an A to D conversion (-32768 < fj < 
32767). If so, there would be 10 sensor parameters per coil (aj , bj , and M 
= 8 cj m ) and 40 sensor parameters to determine. If L = 0; then M = K = 7, 
C j m = wj m , Fj m (uj ) = Uj m and there may be but 36 sensor parameters to 
determine . 
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The measured fluxgate field values (3a) are, within the precision of 
digitized measurement (a j = afj = +0.5 for j i 4), theoretically 
indistinguishable from the result of a system response matrix />jk operating 
on the true field components Bk(x*j) 

3 . 3 

y\ — E />jkBk(x*j) = E Gjk(B,a;h)ajkBk • ( 4a > 

J k=l k=l 


The matrix notation for (4a) is y = pB. The measured scalar field (3b) is, 
within the absolute accuracy of measurement (C 05 from (1)), theoretically 
indistinguishable from the true scalar intensity 

y 5 = [ E (B k (x*5)2]l/2 . (4b) 

k=l 


In (4a-b) 

Bi c (x*J) = field component at jth sensor in bench coordinates. 

ajk = cosfljk and 0 jk is the angle between the jth sensor axis and the 
kth APAFO bench coordinate system axis. 

3 

Note E (aik)2 = l, so 
k=l 


an = 

(1 

- a 12 2 ‘ 

- ai3 2 ) 

1/2 

a 33 = 

(1 

- a 31 2 

- a32^ ) ^ I ^ 

a 22 = 

(1 

- a 21 2 ‘ 

- a23 2 ) 

1/2 

a 43 = 

(1 

- a 41 2 

- a 42^)^^ 


Gjk accounts for possibly non-linear cross-coupling and may depend upon 
the 3 elements Bp of B(x*J), the elements ajk of matrix a, and 
parameter matrix h (see Appendix A). The functional form of 
Gjk(B,a;h) and the numerical values of h must be specified by 
ground calibration. For APAFO we anticipate linear cross-coupling, 

Gjk - 1 and pjk - ajk (M. Acuna, 1990 pers. comm.; see Appendix A). 

is the 4 x 3 system response matrix (see Appendix A). If Gjk depends 
on Bfc (x*J , t i ) , then / 3 jk is implicitly time dependent; this is not 
anticipated, but the possibility is admitted here. 

Now 4 of the 12 ajk would follow from (5) and values for the other 8 ; yet 
all 12 are needed. Thus, 40 + 12 = 52 true calibration parameters need 

to be estimated along with Bk(x*J l ti). (See Appendix A; for linear 
response, there would be 36 + 12 = 48 parameters to estimate.) 

With (3a-b), measured values yj calculated from the bits, the true 
values of the sensor parameters, and C would be true values to within 
instrumental precision; with (4a-b), true values yj follow from true values 
for />jk and Bk(x*J). True values cannot be assumed; yet given (3) and (4), 
the difference between true and measured values can be minimized provided 
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the calibration parameters are steady during each calibration run and 
provided an expression can be identified which links B(x*j,ti) to B(x* 3 ,t£) 
for j £ 4. 

Assume that the ambient field is, at any t, homogeneous across the 
bench. Then the bias field (defined for j i A), 

6 ( x* j , t ) = B(x*5, t ) - B ( x* J , t ) (6) 

is due solely to B? + B 3 . If 6 were steady, then it would be 
indistinguishable from steady VFM biases and no extra steps are needed. To 
see this note B(x* 5 ,tj.) would be B(x*j,ti) + $(x*J), so at each ti 
3 4 

Bfc(x*5) = E (p^-l-kl.E /’jiyj + <5k(x*J) 

1=1 j=l 3 

would be indistinguishable from adding bj* = E aj _ 3 - ^jn^n to b j s 

3 4 n=1 

B k (x*5) = E (pip)- 1 kl E />jl(yj + a j b j*> 

1=1 j=l 

where matrix inversion preceeds assignation of element indicies. 

Now 5 may be steady to within mission goals during certain select 
intervals (preferred VFC times), but might not be at all t. Indeed, 6 will 
not be steady enough when time-varying spacecraft field gradients are too 
steep at the end of the boom. For ESA/PPF torquer rods amounting to a 
point dipole of maximum moment 495 Atm 2 , the radial gradient of the radial 
field over 0.5m amounts to differences of at most, 3. 9 nT at 14 m (the 
radial field is at most, 36 nT; see, e.g., Figure 2). The transverse 
gradient of the radial field, the radial gradient of the transverse field, 
and the transverse gradient of the transverse field, respectively, give at 
most 1/3, 1/2, and 1/6 of this extreme value. In the more realistic case 
of a dipole of moment (350, 0, -350) Atm 2 at the origin, the field at (12, 
7.5, 0) m differs from that at (12, 8 , 0) m by (1.7, 0.5, 0.7) nT, but only 
1.6 nT scalar. 

Corrections for time-varying bias fields of such magnitude are needed. 
The corrected value of fj would be measured fj minus Ek a j”^/*jk^k --with 
all quantities evaluated at ti- Such corrections require values for aj , 
possibly time-varying p, and time-varying 6 in the bench coordinate system. 
The latter require information on the spacecraft field and on the relative 
orientations and positions of the bench and spacecraft coordinate systems— 
perhaps from outboard and onboard star cameras. Now 14 m from a 495“ Atm 2 
dipole-a 0 . 5 ° relative orientation error in a spacecraft field of at most 
36 nT implies a spacecraft field error of at most (36 nT)Xsin0.5° = 0.3 nT ; 
the implied error in the bench bias field correction is negligible (being 
at most ( 3 . 9 nT)Xsin0.5° = 0.034 nT) . Such small rotations through the 
spacecraft field need not alter bench bias corrections appreciably, so it 
is argued that these particular corrections may presume a rigid boom. Here 
it is first supposed that 6 is steady, but that step (II) will allow a 
correction for the part of 6 due to the time-varying primary spacecraft 
field. Corrections for the part of 6 due to bench instrumentation might be 
achieved using ground calibrated models for said instruments. 


10 


Because sensors j = 1 and j = 4 have nearly parallel axes, there are two 
quasi-independent VFM calibration configurations. Configuration 1 uses j = 
1,2,3 (and 5). Configuration 2 uses j = 2,3,4 (and 5). Baseline 
Configuration 0 uses j = 1,2, 3, 4 (and 5). Each configuration can be used 
to estimate the calibration parameters and the uncorrected geophysical 
parameters (BA + B s + B p ) . Comparison of these results enables: 
validation, estimation of calibration uncertainty, and a check on formal 
error estimates. Such comparisons may also prove sensitive to a selected 
(z) component of bench fields due to the 7 -cm separation of coils 1 and 4. 


With the understanding that bj is now the bias due to both the j th 
sensor and any uncorrected steady bias field, and with Bk now regarded 
Bfc(x* 5 ) f the 3 configurations give: 


Configuration 0 
(j = 1.2,3, 4) 

Configuration 1 
(j = 1,2,3) 


yi = E implies Bk = E (^ T ^)" 1 kl > E />jiyj 

J k=l J 1=1 J=1 

3 3 3 

yi = E implies Bk = E ( 1 /> _1 )ki.E Vjiyj 

J k=l 1=1 1=1 


Configuration 2 
(j = 2,3,4) 


3 

y i = E 2 />ikBk implies Bk 
J k=l 


3 4 

= E ( 2 />) - 1 kl E 2 ^jiyj 
1=1 r j=2 


as 


where l/?jl and ^P]l are the relevant 3x3 sub-matricies . Only Baseline 
Configuration 0 is treated here because all calibration parameters are 
obtained and because Configurations 1 and 2 are easier; moreover, sensors 
(1,4) have the strongest cross-coupling and are but quas i- independent . 
Either Configuration 1 or 2 suffice because the resulting Bk at many ti 
overdetermine the calibration parameters for the remaining sensor; yet the 
opportunity for validation offered by multiple configurations is not to be 
overlooked. It seems wise to mount the VFM so that the nearly parallel 
sensors detect mostly the geocentrically radial ambient field component. 


Assume simultaneous measurements at I times ti. With the functional 
forms of Fj m and Gjk specified, and measured bits uj at times ti denoted 
uij from telemetry, we need to estimate 

(i) 52 steady calibration parameters p n = {aj , bj , cj m , p jk) 

(ii) 31 uncorrected field parameters Bik = Bk(x* 5 ,ti) 


from the single scalar calibration coefficient and the 51 measurements: 
(fij ,uij } (41) and fi 5 (II). As sensors j = 1 and j = 4 are nearly 
parallel, it seems there are but 41 pieces of information to estimate the 
31 + 52 parameters. So I > 17 is essential, and I - 10^ is anticipated. 
The 51 equations are 


M 

yij = a j t^ij _ b j ] + E cjm F jm( u ij) (j ^ 4) 

m=l 

M _ 3 

- a j [fij - bj ] + E cj m Fj m (uij) = E />jk B ik = Yij 


(7a) 
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(7b) 


3 

7i5 = Cf i5 = [ E (B ik )2]l/2 = yi5 
k=l 


where the approximations hold only to within the precision (7a) or absolute 
accuracy (7b) of measurement--and only then if both the theory on the 
right and the instrument model on the left are correct. The precision of 
measurement is accommodated by weights. 

For the optionally constrained, weighted, iterated, linearized least- 
squares approach anticipated, we need: measured fij ( j=l, 2, 3 , 4 , 5) and uij , 

C, the functional forms of Fj m and Gjk; the uncertainty estimates o\ in the 
measured fij (again, nominally +0.5 bit for j < 4, and a5 from (1)), and 
some initial guess of the calibration parameters--the measured values from 
the ground calibration. The idea is that minimizing the difference, albeit 
weighted-sum squared, between measured and theoretical true values of yij 
estimates true yij . 

Following Lancaster et al., focus on f ij rather than yij* From (7), the 
theoretical true values fij are 

3 M 

fij = E ^jk B ik - E c jm I? jm( u ij ) ] + bj (j ^ 4) (8a) 

k=l m=l 

3 

fj.5 = C-l[ E (B ik ) 2 ]l/2 . (8b ) 

k=l 


The partial derivatives of fij with respect to calibration parameters p n 
and uncorrected geophysical parameters Bik are, respectively 


3 f ij __ 

= s ijn 

9pn 


0f 


1J 


9B ik 


= r ijk 


with subscripts i (1+1) for time, sensor j (1+5), calibration parameter n 
(1+52) and field component k (1+3) in bench coordinates. 


Now 9fi5/0Pn = ^i5n» so the only non-trivial Sij n have j i 4 and, by 
(8), are 


9fij 3 M 

= £ /?jk&ik - E c jm F jm( u ij)] = ~(fij - b j ) / a j 

3aj k=l m=l 


= 1 

9bj 

9f ij 

= -aj-lFj m (uj.j ) 

9 c jm 


(9a) 


(9b) 


(9c) 
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(9d) 


ij 3 8/>j 1 3 

= aj-1 E [ Bn] = E aj“l<5kl B il - a j “ . 

9/>jk 1=1 9/>jk 1=1 

Note (9d) anticipates linear cross -coupling and treats the elements of p as 
independent. For non-linear cross-coupling (9d) would be replaced by 

9fij 3 9/?j i 3 9a ji 9Gji 

= aj -1 E [ Bn ] = E aj-l{[Gji( ) + ( )aji]Bii> 

9aj k 1 =1 9aj k 9ajk 9ajk 

where some of the 9aji/9ajk are zero by (5). For no cross-coupling, all 
Gji would be 1 and all 9Gji/9ajk would be zero; for non-linear MAGSAT 
cross-coupling ,3 of 9 Gj i differed from unity, 9Gji/9aj’k was zero unless j 
= j’ — in which case 9/Jjl/9ajk depended only upon Bi and the ground based 
values of h. The replacement for (9d) would thus suppose j = j’ and 
estimate the 8 independent ajk rather than 12 />jk* 


The only non-trivial rijk are, according to (8) 

9fij 3 9pjl 

= aj-1 E [/>ji<5ik + ] Bi i (10a) 

9Bik 9Bik 

9fi5 

= (C2f i5 )-lB ik . (10b) 

9B ik 

The partials (9-10) imply sij n = sij n (p n , Bii ; f ij , Fj m (uij ) ) and rijk = 
rijk(Pn» B ii; f ij • F jm(uij ) ) • The inverse problem is thus non-linear: refined 
estimates of p n and Bik require prior values for p n and Bik» In the 
iterative linearized method of estimation anticipated, all partial 
derivatives are evaluated at the best previous values of p n and Bik (and 
f ^ j ) — be it an initial guess or a subsequent estimate. 

For quantity q, denote: true value q; guess q; measurement q; and 
estimate q. The expected VFM change upon launch enjoys zero mean, so 
initial guesses of aj , bj , cj m , and aj k are pre-launch values ’measured* 
during ground calibration (p n = 'pn) . For purely linear cross-coupling an 
initial guess of Bik follows from p n and telemetry. To obtain Bik for non- 
linear cross-coupling, guess B°ik (e.g., suppose either linear response or 
the values indicated by position and attitude in a 29,500-nT geocentric 
axial dipole) and use it with ajk and h to compute G°jk, hence />°jk- Then 
use /j°jk and yij to compute B^ik and iterate until satisfactory values for 
Bik and Pjk are obtained. Either way, these guesses give 


yij 


#v — rw M n — 

a j [ f i j - bj ] + E c jm F jm( u ij) 
m=l 


(j * A) 


(11a) 


/ 
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(lib) 


3 ~ ~ A ~ 

Bik = E [^ T / , ]" 1 kl E />jiyij 
1=1 j=l J J 

3 - 

f i5 = C-l[ l (B ik )2]l/2 f f i5 (He) 

k=l 

where the anticipated inequality on the far right of ( 11 c) is to be 
eliminated, along with the guesswork, by calibration. 

Let f = { f i j w fi 5 > be the column vector of 51 true fine readings, z = 

{p n , Bik) the vector of 31+52 true calibration and total field components, 
and let Y = {sij n >rijk} t> e th e matrix of the partial derivatives 9f/9z 
(evaluated at true z) with 51+52 rows and 31+52 columns. Let V be the 
51x51 matrix of fine measurement weights--the expectation value {G'ijfl’ik )" 1 
(noting that the Jjj are expected to be independent and steady so that W is 
a diagonal matrix of constant element^) . 

We seek to minimize (f - f) T W(f - f ) . Unfortunately, constraints may be 
needed to hold estimates of a few p n near previous estimates (if not ground- 
based values). The reason is that some data bits comprising some uj are 
not expected to change independently during some calibration runs due to 
VFM attitude in the largely dipolar ambient field. This problem is 
exacerbated by the loss of scalar signal (LOSS) encountered when the field 
vector is inclined by more than 60° relative to the single SHM cell axis. 

Two suitably oriented SHM cells will prohibit LOSS events during intervals 
of nominal spacecraft attitude. During such runs, the corresponding p n (a 
few of the 28 wjk) cannot be estimated. Either these special parameters 
must be fixed to their previous values and the corresponding rows and 
columns of T deleted or, perhaps more simply, a constraint can be added to 
the objective function which keeps them near their previous values. The 
added constraint is of the form (z - z)A(z - z) where z is the previous 
est imate--regarded as a ’measured’ value. "Nearness" is judged by the 
elements of the damping matrix A which, in turn, is a modification of the 
information matrix from the previous calibration 0 “ 1 . 

The previous calibration covariance Q is generally a full matrix and is 
the expectation value {(<7z)(<7z) T }, but only its upper left-hand portion 
contains the expectation value {(ap)(tfp)T) (i.e., the covariance of p from 
the previous estimate). The information matrix Q - 1 will normally be the 
expected weight or confidence matrix. However, if the parameters change 
too much, this weighting may need to be modified or reduced. In the 
extreme case of drastic change from the previous calibration one would 
assign zero confidence to all previous estimates of B^k and Pn --excepting 
those which cannot be estimated due to LOSS or attitude restrictions and 
are, perforce, kept close to the previous estimate. This can be achieved 
by setting all elements of Q-i to zero except the few elements in the upper 
left corresponding to the targeted (^Pn^Pm)* For example, if only P 7 is to 
be constrained, then only A 77 = X{ap 7 ap 7 >“l would be non-trivial. 
Alternately, if all p n are to be constrained, then the non-trivial elements 
would fill the upper left-hand portion of A. The damping parameter X < 1 
may reflect reduced confidence in the relevance of the previous 
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calibration. We aim to set X » 0 during the initial, on-orbit check-out 
and on other selected occasions. It is also possible that we may instead 
use the fact that some p n enjoy hard prior bounds (e.g., I a j k I ^ 1 ) • 

In view of ( 2 ), a set of ti which fails to span zero-crossings of the 
offset field (the 7 -bit uj changing from 1111111 to 0000000) will not 
sample a flip of at least one ujk (i.e., Uj7); then at least w j 7 is not 
resolved. With appropriate orientation of the SHM axis, resolution might 
be achieved by carefully selecting the t£. If a polar orbiting sensor (say 
j = x) is always nearly parallel to the geomagnetic X component of a 
roughly axial dipolar field with a -30000-nT equatorial field, then both 
u x 7 and u x g will flip — but they will always flip together. Then 
calibration coefficients (w X 7, w x g) are not separable as they cannot be 
determined independently. This indeterminacy might be considered 
irrelevant: a single determinable calibration coefficient could replace the 
two. A similar situation exists for Wy7, Wyg, perhaps Wy5, and conceivably 
Wy4; no such situation is expected for Z. 

Determination of the set of all calibration parameters requires exercise 
of all relevant offset field generator bits for all sensors; unfortunately, 
it may not occur at times when scalar data is available. Full exercise of 
all bits with scalar data is requested periodically to separate and 
identify drift in all wjk* More frequent exercise is desirable to 

determine more nearly current values of Q. If the x,y t z sensor axes are 

nearly parallel to the X,Y,Z geomagnetic components (as are the nominal 
spacecraft coordinates when ascending above the Northern Hemisphere ) , then 
full exercise is expected to occur naturally on the z (j=l,4) axis, might 
not occur naturally on the x (j=2) axis, and will not occur naturally on 
the y (j=3) axis. As to LOSS, consider a polar orbit through an axial 
dipole field. If the SHM axis points in the x (ram) direction, the LOSS 
cones would be poleward of +40.9° latitude; if the SHM points in the z 
(down) direction, the LOSS band in such a field would be equatorward of 
+16.1°; if the SHM is oriented at 45°, the LOSS zones in such a field would 

span 61.8°N to 7.6°N in the descending (day) arc and 61.8° S to 7.6°S in 

the ascending (night) arc. The 98°, rather than 90°, orbital inclination, 
11° dipole tilt, and non-dipole field will shift the LOSS zones. The first 
effect bars data poleward of +82° --suggesting an off-nadir SHM attitude; 
the second effect would seem to narrow any equatorial LOSS zone. A 
detailed simulation seems appropriate in view of non-dipole fields. 
Depending upon the SHM orientation, some separation of the more significant 
wjk may be achieved by attitude maneuvers: a 90° roll at high latitude 
likely results in LOSS, but a 45° roll may allow separation of many Wyk . 

With z denoting the previous value for z, the total objective function 
to be minimized is 

A 2 = (f - f)T W (f - f) + (z - z)A(z - z) . (12) 

The linearizations 

z - z (12a) 

f 2 f = f + T(z - z) (13b) 
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yield f- f-f- f=f- f + T(z - z) and 

~ -A " ~ "A 

&2 = [(f - f) + T(z - z)]Tw((f - f) + r(r - z)] 

+ [(z - z) - (z - z)] T A[(z - z) - (z - z) ] . 

A 2 a s approximated by (13c) is minimal if, and only if, 

r T W[(f - f) + T(z - z)] + A[(z - z) - (z - z) ] = 0 


(13c) 


(14a) 


or, with d = f - f 

z = z + [Pwr + A]-l[nwd + A(z - z)]. (14b) 

It is stressed that the vast majority of the elements of A are zero — 
notably all those corresponding to any value for Bifc. The few non-trivial 
elements of A corresponding to the unresolved p n also correspond to guesses 
z which are equal to the previous value z. So A(z - z) = 0. The inversion 
of the 31 + 52 square matrix [f T Wf + A] is formidable with A = 0, but must 
be done in occasional full-range calibrations. The routine inversion with 
non-trivial upper left A from the pxp sub-matrix of Q _ 1 is described by 
Lancaster et al.,1980; then inversion of even the 31x31 sub-matrix 
associated with the Bik is avoided by partitioning it into I, 3x3 
matricies . 

With the first estimate of true z, we can estimate T = 3f/9z|~ , 

~ I* 

set d = f - f, and calculate the refined estimate 


l = i + [Hwr + A*]-i[Pwd + A*(z - Z)] . ( 15 ) 

Note A* may differ from A because Q*”l may differ from Q~1 in view of z and 
the likelihood that exposure to the space environment cause p n to deviate 
from previous values. Indeed, z itself in (15) might be replaced with z or 
higher iterates if they prove that p n evolves due to fluxgate drift and 
cosmic radiation. Of course , weights W are dictated by instrumental 
precision and are thus fixed. Extreme care is required in the selection of 
A*. If estimate z from (14b) differs enough from guess z, then estimate f 
can be far from guess f; similarly, if estimate p n differs enough from 
previous value p n , then the previous Q can be a poor estimate of current 
calibration covariance and both p n and A, worthless. The former case 
demands iteration; the latter case indicates attempting to exercise all 
bits, setting Q*-l and A* to zero, and redetermining all p n . Clearly, no 
confidence in prior estimates of the calibration parameters is warranted if 
a micrometeor impact reorients the instruments! 
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Phase One: Step II, Spacecraft Field Parameterization 

This section outlines a parameterization for the spacecraft fields, both 
primary and secondary. This could readily be translated into a least- 
squares problem. Included in the parameterization are the torquer rod 
fields, fields from the solar array, and fields induced both in conductors 
via Faraday’s law and in magnetically permeable materials by ambient 
fields. Further study will be undertaken in Phase C to refine this 
parameterization and to translate it into operational software. 

However, the outlined procedure may be unnecessarily complex. For 
example, it is hoped that induced fields will be negligible and that the 
solar panels will be so wired as to eliminate their magnetic field. In any 
case, a model including the magnetic torquers and, probably, the solar 
panels will be required. For the present discussion, a model formulation 
for induced fields will be included but not practically developed. To 
determine if they may be ignored during the mission, an appropriate 
magnetic survey of the spacecraft should be conducted prior to launch. 
Otherwise, a model such as the one presented must be implemented prior to 
launch to include this contingency in our plans. 

In addition to these assumptions are the following: 

1) Torquer rod, solar panel and other currents will be available in 
the telemetry with sufficient accuracy and time resolution to permit 
spacecraft field calibration and application of that calibration to 
correct the measured field. 

2) Changing states of spacecraft and instrument systems will be known 
via telemetry. 

3) Boom motion can be ignored. As will be shown, this assumption is 
justified as long as boom angles relative to nominal are less than 
0.6°, and possibly if they are less than 2.5°. Boom model studies 
should be conducted to confirm that these limits are not likely to be 
exceeded. Note that if boom motions do exceed these values, but do 
so in a way that can be geometrically modeled ( i . e ., following the 1st 
normal mode), then the combination of stellar attitude measurements 
on the main spacecraft and on the APAFO optical bench should allow a 
correction for boom motion. 

4) It is assumed that the torquer rod field strength depends linearly 
upon the applied current with no appreciable hysteresis. This 
implies the use of ultra-soft, low-hysteresis rods. If necessary, 
the model could be made to account for small non-linearities. The 
rods must be checked pre-launch for linearity and hysteresis. 

Given these conditions, it should be possible to calibrate the fields 
from the torquer rods and solar panels to within 12 and the entire 
spacecraft field to within 52. 

Many spacecraft field paramete rizat ions are possible, but the purpose of 
the spacecraft field model is to enable accurate calculation of the ambient 
field by correcting VFM calibration and the total field measured. This 
implies that the spacecraft field parameterization must represent the 
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spatio-temporal variation of real spacecraft field sources. We anticipate 
a full complement of spacecraft field sources: bi-axial torque rods rated 
at 350 Atm 2 each; solar array panel current; other designed currents; 
currents induced in the spacecraft and the APAFO equipment (boom, bench, 
etc.) by the time-varying spacecraft and ambient fields; induced 
magnetization; and permanent magnetization. The torque rods are considered 
primary sources; the solar array and induced fields are considered 
secondary sources. Figure 2 shows the spacecraft field as a function of 
boom length for two torquers of 350 Atm 2 in the configuration of Figure 1. 
Other sources may be weak, but must be anticipated unless pre-launch 
measurements of the field near the spacecraft indicate they are indeed 
absent. Ground-based spacecraft field calibration must determine the 
relations converting torquer rod and solar array current telemetry to the 
magnetic field produced at the expected bench location. 

Application of Step (I) yields initial estimates of the total magnetic 
field at the outboard location in bench coordinates at times t^. It also 
gives calibration parameters for estimating the field at other observation 
times. At any observation time tj the spacecraft field in bench 
coordinates is 

S3 S 

B jk = E *jkl B ’jl (16) 


where B^’j! is the 1th component of the spacecraft field at tj in 
spacecraft coordinates and the Rjkl are elements of the matrix R(t n ) which 
rotates the spacecraft coordinate system into the bench coordinate system 
at t n . Clearly R is time dependent due to boom motion. The reverse 
rotation of bench coordinates into spacecraft coordinates is denoted S(tj). 
The Rjkl can b e derived from the attitudes determined by the bench star 
cameras and by the ESA/PPF onboard star cameras. ESA/PPF attitude 
knowledge is put at 0.01° (one sigma), but pointing accuracy at the torquer 
rods is put at 0.033° during disturbed times (120 arcseconds, or about +0.6 
mm positioning at the ends of a 1 m rod); the star camera axis to VFM 
optical cube pointing accuracy is put at better than 12.5 arcseconds rss 
(perhaps 5 seconds). So the relative orientation of the torque rods and 
bench coordinates should be determinable to about 121 arc-seconds rss. 

At present, it seems that boom motion through the spacecraft field (as 
distinct from the ambient field) might be a negligible source of error due 
to the small amplitude of the spacecraft field at the bench. The rigid 
boom approximation, 

S 3 S 

B nk - E Rk 1 B * n 1 (17) 

1=1 

where the R^l are the mean values of the Rjkl» ma Y thus prove useful. 

To see this, recall that for the realistic example in section (I), a 
"torque rods" magnetic dipole of moment (350,0,-350) Atm 2 at the origin 
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generates a field at (12,7.5,0) m which differs from that at (12,8,0) m by 
(1.73,0.50,0.69) nT, but only 1.63 nT scalar. The field at (12,8,0) is 
(12.56,16.15, 11.67) nT, or 23.56 nT scalar. These realistic values 
suggest that the 36-nT worst-case spacecraft field estimate is a little too 
pessimistic. Also recall that a 0.5° relative orientation error in a 36-nT 
spacecraft field gives at most a 0.3-nT error in the spacecraft field 
estimate. This approaches the 0.5-nT APAFO goal. Now, relative orientation 
fluctuations of +0.5° due to flex of a parabolic 12-m boom imply +5 -cm 
translations of the bench through the spacecraft field. These in turn 
imply a spacecraft field correction about +10Z of those noted above for the 
0.5-m displacement, or about 0.4 nT in the worst case. Because errors of 
0.3 nT and 0.4 nT yield 0.5-nT rss, fluctuations in the relative bench- 
spacecraft orientation of up to 0.5° are clearly magnetically tolerable 
(i.e., negligible). In contrast, fluctuations over + 1° may prevent 
meeting the 1.0-nT spacecraft field requirement. 

In fact, a 2.86° fluctuation in the relative orientation (known to +121 in) 
of the bench to the spacecraft due to boom flex would imply a 30- cm 
translation at the tip of a parabolic 12-m boom (60 cm for a linear boom) . 

In the realistic example, this would imply a 1.0-nT correction to the x 
component of the spacecraft field. But the parabolic boom model cannot be 
assumed reliable because the orientation fluctuation may be caused by an 
overtone rather than the fundamental, roughly parabolic, flexural mode. 

Thus the desired correction cannot be made with certainty and the 1.0-nT 
spacecraft field requirement cannot be met. Clearly if elements of Rjkl 
differ from the mean Rk]_ so much as to imply a +2.8° fluctuation, the data 
must be flagged for a potentially serious, off-nominal relative bench- 
spacecraft orientation. The point here is that even +0.5° fluctuations 
seem large, yet are quite tolerable; thus the rigid boom approximation may 
in fact prove quite useful. 

Parameterization. The chief spacecraft field at a position x in 
spacecraft coordinates is here represented as 


rods solar in 

B s, jl(x) = Bji(x) + Bjl(x) + Bji(x) + B 0 i ( x) (18) 

where Bji is the induced field and B 0 1 is the steady spacecraft field due 
to both permanent magnetization and steady currents. 


Torque Rods . In spacecraft coordinates the magnetic field in nT at x 
(in meters) generated by the thin solenoid (two point monopoles) model of a 


torque rod with fixed ends at x+ and x~ s 


eparated by D = I x* - x~ I 


is 


rod 

B (x) 


10 2 NIA 


D 


x - x+ 

[ 



(19a) 


where N is the number of turns, I is the solenoid current in Amperes, and A 
is the cross-sectional area of the rod in m 2 . NIA is the absolute Amperian 
magnetic dipole moment. The field due to two solenoids is, at time t n , 
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+ 

rods 2 10 2 N r A r xi - x r i xi - x r i 

B nl (x) = E [ - ] l nr ( 19 b) 

r=1 D r lx - x +|3 |x - x -|3 

2 r r 
= E M lr (x)I nr . 
r=l 

If x is time dependent, then so is Mi r (x); such dependency would be omitted 
in the rigid boom approximation. Prior values for Mx r at x* may be taken 
from ground-based design and pre-launch testing, but estimates of the 6 
elements of H need to be made in-flight. 

The relations between the true magnetic dipole moments and the torque 
rod currents are not known and must be specified; at the very least, a 
soft, high -permeability rod core implies that the true magnetic dipole 
moment will greatly exceed the vacuum core value unless either N r A r is 
regarded as an effective turn-area or I nr includes the effect of 
magnetization. The linear relation 

inr - a r(inr - Pr) ( 19 c) 

may prove adequate for elect romagnet- type rods--with ar regarded as the » 
permeability gain and /? r the offset or bias of rod r; a cubic or higher 
order polynomial will be used as needed. Hysteresis in the rod cores may 
be observable, but no correction for it is planned due to the complicated 
time history of rod currents. It seems that ultra-soft, low-hysteresis 
rod cores are required; even then, a hyperbolic tangent may be needed to 
mimic an approach to rod core saturation. 

Solar Array . In practice, the solar array fields could be made 
negligible by suitable back-wiring. If not, then a model is required. For 
a rotating solar array, the angle of rotation will need to be known for the 
model. If instead of complete back-wiring, adjacent panels in the solar 
array are intended to cancel one another, additional information will be 
needed. In particular, the multi-panel geometry should be specified and 
the telemetry signal should give information about which panels are 
operable and, if possible, what current is output from each panel. 

The solar array wiring and current flow design are not yet known, so 
detailed modeling of the field produced by the solar array is not yet 
possible; however, we anticipate a representation in terms of two or more 
rectangular current loops. For example, in the special case where the 
solar array field reduces to the 2 counter-circulating current loops of' 
area A4 which are far ( I X-X4 I 2 » A4) from the bench, the array field is 
well approximated by two anti-parallel point magnetic dipoles of moments 
+KIj4A4m located at xi and X2 : 

solar (ni*m)ni - m (n2*m)n2 - m 

B(tj ) = 102KA A ( - ]I i4 (20a) 

I x - xi I 3 I x - X2 I 3 
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where ni = (x-xi) / lx-xi I and n2 = (x-X2) / Ix-X2 I are the unit-normal vector 
point from the dipoles towards the field point at x, and where K is the 
fraction of I j 4 flowing around the loops. Now m, ni, H 2 t xi, and X2 are 
all time dependent because the solar array rotates relative to the 
spacecraft so as to face sunward; yet ’estimates’ for all these quantities 
can be derived from ESA/PPF design, orbital data, and ephemeris data. 

More generally, if a linear relation between true and telemetered solar 
array current is presumed, then, at time t j , 

solar _ 

B j 1 (x) = Q j 1 1 j 4 ® Qj 1^4 ( j 4 • P 4) • ( 20 b) 

Nominal Q j 1 is expected to have a distinct, ephemeris specified time 
variation which repeats every orbit; special calculation of the Q j 1 is 
needed for off-nominal attitude. Refinement of these calculated Q j 1 is not 
anticipated; instead we expand the coefficient <*4 into the 3 component <*14: 

solar _ 

Bji (x) = QjiaiAdjA - PO • ( 20 c) 

In-flight estimates of ai4 which deviate significantly from a single 04 
indicates inadequacy of either the solar array model leading to the form of 
Q j 1 or in the prior estimate of it. 

Induced Fields . In the spacecraft coordinate system, Faraday ’ s law 
implies that the time rate of change of ambient, torque rod, and solar 
array fields are equivalent to the negative curl of an induced electric 
field. This induced electric field drives electrical currents in every 
electrical conductor onboard the spacecraft. These currents generate the 
"Faraday" induced spacecraft field which, by Lenz’s law, opposes the 
inducing field. In addition, the spacecraft and ambient fields induce 
magnetization in permeable materials onboard the spacecraft (notably the 
torque rod cores); this induced magnetization generates the "Henry" induced 
spacecraft field. 

Ideally, an explicit model of Faraday induced fields might begin with 
the position and electrical conductivity a(x) of every conductor onboard 
the spacecraft and an estimate of the electric field E induced by the total 
space- time-varying magnetic field; proceed through evaluation of the 
current density J induced in the conductors by the induced electric field 
and calculation of the induced magnetic field via the Biot-Savart law; and 
wind up in iterative correction of the total time-varying field for the 
Faraday induced field. Similarly, an explicit Henry induced field model 
might begin with the position and permeability Jl(x) of all magnetizable 
materials onboard the spacecraft. Alternately, one might solve the vector 
magnetic diffusion equation (e.g., 9 tB = lB for linear J = aE and 

B = /IH). The short magnetic diffusion times expected for thin ESA/PPF 
conductors indicate a Faraday induced field which is in phase with the 
time-rate of change of the inducing magnetic field at the ultra-low 
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frequencies of interest. Similarly, magnetically soft materials will have 
an induced magnetization, and Henry induced field, which is in phase with 
the ULF inducing magnetic field. 

Only the induced field at the bench is of interest, so we bypass the 
exhaustive explicit procedure by supposing that the Faraday induced field 
at x consists of a part proportional to the time rates of change of: the 
torque rod field (which is in turn proportional to the torque rod current); 
the solar array field; and the ambient field. Similarly, the Henry induced 
field may be bypassed in favor of a simple transfer function. In this 
"lumped element" approach to induced spacecraft fields, the 
parameterization of the induced field at the bench is 

in 2 A 

Bji = E Di r d t Ijr + Elk^tBjk + Fi4[Pji9 t Ij4 + (OtPjlUjAl 
r=l 

2 3 A 

+ E Ti r Ij r + E Ui k B jk + V 14 CjiIj 4 (21) 

r=l k-1 


here 9tljr * (lj*r - Ij-ltrJ/Ctj ~ * 0 - 1 ) ~ fl r®tljr/9t anc * 9t^j4 = 

(Ij, 4 - Ij-l» 4 )/(tj - tj_i) = fl 4 9lj 4 /9t. The ambient field in the 
spacecraft coordinate system can initially be taken as the rotated total 
field from the calibrated VFM: B A j k - EnSjmnBjnCx^) ; similarly, 9tB A j k - E n 
S*jmn(Bj , n (x 5 ) - Bj , n (x 5 ) ] / ( t j - tj_i). (Although the rigid boom 
approximation may be adequate, use of time -dependent S(tj) reduces 
contributions from boom twist through the ambient field; corrections to 
9 t B A j k due to the time-varying spacecraft field at x^ can be included on 
higher iterations). In principle, the elements of (D, E, F) could be derived 
from the Biot-Savart law by integrating over the result of geometric 
operations on the induced current density; similarly (T,U,V) would follow 
from the curl of a dummy integral over the result of a geometric operator 
on the induced magnetization. P j i and C j i account for the time variation 
of the solar array orientation. For simplicity, Pj i is taken as Qj \ 
evaluated at an x near the base of the solar array boom (where the solar 
panel field is much stronger than at the center of conductivity); C j i is 
taken as Q j i evaluated at an x near the torque rods (where virtually all 
the permeability is located). In the unlikely event that a spatial 
parameterization of the induced field is needed for bench gradient 
corrections, fields mapped by (D , E , F) might be treated as follows: (D) a 

point dipole located closer to the torque rods than to the center of 
conductivity of the spacecraft; (E) a point dipole at the center of 
conductivity; (F) a point quadrupole much closer to the base of the solar 
array boom than to the center of conductivity. The fields mapped by 
( T , U , V ) might be treated as: (T) two-point dipoles near the center of 

permeability; (U) a point dipole at the center of permeability; (V) a point 
quadrupole between the solar array boom base and the center of 
permeability. It might further be supposed that the center of permeability 
would coincide with the torque rods, while the center of conductivity would 
coincide with the center of mass. 
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The spacecraft field model parameterization is now 


S 2 _ - . 

B ’u = E Mi r a r (Ij r - p r ) + Qjl<*14<Ij4 - PO 
J r=l 


E Di r a r 8 t Ijr + E Eik9 t Bjk + *14 t p j l“140tlj4 + (8t*jl)«14(Ij4-/>4>] 
k=l 


r=l 

2 


+ E Ti r a r (Ii r -^r) + E UikBjk + V^Cj iau (Ij 4"^4) + B ol • 
r=l k-1 


(22a) 


The effects of Mi r are virtually indistinguishable from those of Ti r , as 
are those of Qjl> Fi A 9t p jl> and ^ 1 AC j 1 1 

B S, jl = E [Mi r +Ti r ]a r (Ijr - Pt) + tQj 1 +F 14 (9t p j l> +v 14Cj i ]ai4 ( I j 4 - PO 
r=l 

2 _ 3 A — 3 A 

+ E Di r a r a t Ijr + E Eik0 t Bjk + Fi4Pjiai49 t Ij4 + E U lk B jk + ®ol (22b) 

r=l k=l 


Torque rod offsets p v are indistinguishable from any other steady 
spacecraft field. So, grouping together time-varying rod, solar array, and 
ambient field effects and steady field effects, 

S 2 - — 

B ’u = E {[Mi r + Ti r ]a r (Ij r - P r) + Di r O r (9 t I j r ) > 
r=l 


+ UQjl + Fl4(9tPjl) + v 14 c j 1 ] a 1 4 < Ij 4 - PO + F 14 p jl a 14(9tlj4) > 
+ E (Uik(Bjk) + EikOtBjk)) + (BqI - [Mi r + Ti r ]ar^r> (22c) 


or, more simply 
S 2 

B ’ j i = E {M’ir(Ijr) + D* lr (9tlj r ) > 
r=l 

+ { [Q j 1 + FuOtPjl) + VuCjUaudj, - Po + Fi A Pjiai4(9tij4)> 
+ E {Uik(B^k) + Elk(9 t Bjk)} + { B ’ ol } • (22d) 

k=l 


Here Qn, Pji, 9t p jl* and c j 1 are considered known from solar array design, 
orbital , * and ephemeris data; the command current and its time rate of 
change are known from telemetry; and the guesses for B A j k and 9tB A jk follow 
from step I, VFM telemetry, and the dual star-camera specification of 
S(tj). The 24 spacecraft field calibration parameters TBD are: 6 M’i r for 
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the direct torque rod field and its Henry induced counterpart; 6 D’i r for 
the Faraday field induced by the torque rod field; 3 ai 4 for direct, 3 F 14 
for Faraday induced, and 3 V 14 for Henry induced solar array fields; 3 Uifc 
for Henry induced and 3 E]^ for Faraday induced fields associated with 
spacecraft flight through the ambient field; and 3 components of the steady 
spacecraft field B' 0 i. Note that a suitably wired solar array would 
eliminate 9 parameters and costly evaluation of Q, P, and C. 

Special Calibration Periods . To ease detection of a roughly 30-nT 
spacecraft field superimposed upon a ambient field exceeding 30,000 nT, it 
is useful to subtract an initial ambient field model and correlate 
variations in the residuals with those in Ij r , Ij4, and their time 
derivatives. Unfortunately, the torque rod currents are expected to be 
strongly correlated with the ambient field exerting the designed torque. 

So, as a check on the routine calibration* special calibration periods 
would be desirable for this primary spacecraft field source during which 
the torque rods are pulsed in a preset pattern. 

The pulsing pattern considered is a series of abrupt changes in the 
command current from one level to another level, with the command held 
fixed at each level for an interval At. Start with both rods off. Raise 
one rod to level I Q for At, lower it to zero for At, then lower it to -I Q , 
and finally raise it to zero for At. Then repeat the sequence at levels 
+2I 0 * +3I 0 » ••• ±Imax- Then repeat this pattern for the other bar. With 
At = 1 s and the VFM in campaign mode, induced transients can be monitored 
with 128 VFM samples and at least two scalar samples will be acquired. 
Longer At is needed if transients are a problem. Clearly AT should be much 
larger than the inductive time constant of the torque rods. If the rods 
are pulsed in shadow, the times of zero -level rod current could be ideally 
suited to VFM calibration. 

In the rigid boom approximation, the zero-current VFM measurements would 
be fitted with a low-order polynomial in time; this model would be 
subtracted from the entire measurement sequence, and the 6 elements of 
[Mi r +Ti r ]a r = M* i r estimated by a leas t- squares fit to the residual field - 
which is simply modeled as M’i r Ij r . Rod core remanance should appear as 
residuals at the zero-current times which share the sign of the previous 
command current. 

Boom flex, twist and stretch will contaminate the residuals to a low 
order polynomial due to VFM rotation through the ambient field; this signal 
may well exceed that of the torque rods. If the spacecraft /torque rods 
attitude is known only to 121 " rss, the total field ought not be rotated 
into spacecraft coordinates prior to elimination of some ambient field 
model, lest attitude errors of up to 50 , 000sin(120" ) ~ 30 nT be introduced. 
So, for a flexible boom, first rotate the VFM samples into celestial 
coordinates using the bench star-camera data only ( 12 . 5 " rss); then fit a 
low -order polynomial through the zero -current samples; next rotate the 
residuals into spacecraft coordinates; and finally estimate M’i r from the 
residuals in spacecraft coordinates using least squares. Any prior 
information on Mi r must be distinguished from that on M’i r . 
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The low-order temporal polynomial fitted to celestially oriented, zero 
rod current level, ’ambient' field values might be abandoned in favor of 
simply averaging or otherwise interpolating the (celestially oriented) 
field values during the zero current level times adjacent to the pulse of 
interest. This reduces presumption regarding ambient field smoothness and, 
perhaps more importantly, allows the pulsing to occur at high latitudes - 
where the crossings of auroral current sheets would be poorly represented 
by a low-order polynomial. Pulsing the rods at high latitudes could also 
provide useful zero current level times for VFC calibration. Furthermore, 
hysteretic remanance in the rod cores (which may appear in the residuals to 
a low-order polynomial) might be reduced by suitably weighting the average 
(or other interpolate). 

Clearly 4At must not equal the fundamental period of the boom; selecting 
At equal to this period may improve spacecraft field calibration. 
Alternately, sufficient time should be allowed between torquer steps for 
induced boom transients to die out. Special calibration periods conducted 
with the boom either in constant shadow or constant sunlight should enjoy 
minimal thermal boom transients. 


Phase One: Step III, Iteration 


The initial spacecraft field estimate from (II) is just the torque rod 
direct and Henry induced field B’ s ji ■ E r M’i r Ij r ; in bench coordinates 
this is B s jk = Elr SjklM’i r Ij r , or about Ei r SklM’i r Ij r in the rigid boom 
approximation. To correct the VFM calibration for the time-varying bench 
bias field caused by the torque rods, not only the orientation, but the 
positions xJ = x° + Rx*J of the VFM and the SHM must be known in spacecraft 
coordinates. It can be assumed that shifts in VFM and SHM positions on the 
bench itself would appear as off-nominal calibration parameters or LOSS 
events; indeed, if steady, they might be deduced and used to obtain 
corrected x*J . The key question thus regards the position of bench itself 
x° : did the boom deploy as designed? Clearly, if the Rkl differ 
significantly from designed values, the designed boom twist or deployment 
angle was not achieved - suggesting a deviation from the designed bench 
position. Moreover, if the elements of M*i r differ from the values of 
Mir^r anticipated at the designed bench location by more than explicable in 
terms of Ti r a r , then one may infer that the value of x° differs from the 
designed value. Then the values of the 3 components of x (or rather x 5 ) 
can be calculated from the 6 estimates M’i r using (19b) on the suppositions 
that M’i r = Mi r a r and that both the a t and the torque rod positions are as 
designed. Of course, Rik, and x * 5 are regarded as known, so the 
determination of x 5 implies x° . Thus the position of the bench in the 
primary spacecraft field can be estimated in the unlikely event of an off- 
nominal boom deployment. 
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With M'i r Ij r having the x dependence indicated by (19b,c), and supposing 
Pr - 0 r> the primary spacecraft field can be evaluated at the positions of 
the VFM as well as the nominal SHM position. In bench coordinates, the 
difference between the two values estimates the time varying part of the 
bench bias field 6 . As described in section I, this implies a correction 
to the f i j ; the corrected ffj can then be used to derive corrected 
calibration parameters p n and corrected total as in section I. The 

value of B s j i is then subtracted from the derived using the corrected 

calibration parameters to obtain the next estimate of the ambient field, 
B^jk» in bench coordinates. 
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PHASE TWO 


During phase II the coordinate transformation from the VFM to an earth 
centered system will be confirmed, previously undetected spacecraft fields 
will be isolated, and initial models of the geomagnetic field computed. 

This phase uses procedures developed for reduction of Magsat data. It 
depends upon the fact that the platform is earth oriented and therefore any 
unmodeled spacecraft fields are rotating relative to the field of the 
earth* The basic assumptions are that the magnetometer coordinate system 
does not change relative to the star sensors and that any unmodeled 
spacecraft fields are constant during the period of analysis. 

This phase is accomplished using a spacecraft vector magnetometer 
measurement model in the software used to derive models of the geomagnetic 
field, i.e.,the computer code FIT. 


The magnetic field vector in earth fixed cartesian coordinates is 
represented as 


m 


Bg = Tcs lYl 

(23) 

IzJ 


where X = -B$cos5 - B r sin5 

(24a) 

Y = B^ 

(24b) 

Z = Bgsintf - B r cos5 

(24c) 

(Tcs ) 11 = -cos ( 0-5) cos^; (Tcs ) 22 = cos^; (Tcs)33 = -cos(0-5) 
( T CS ) 12 = -sin^; (Tcs)l 3 = - s in ( 9 - 6 ) cos^ 

(25) 


( t CS>21 = -cos(0-<5)sin0; (Tqs)23 = -sin(0-5) sinfJ 

( t CS)31 = sin(0-5); (Tcs)32 = 0 

B = -VV( r , 0, g n m ,h n m ) , (26) 

where r, 0, $ are the usual geocentric coordinates, <5 is now the difference 
between geodetic and geocentric latitude, and the g n m and h n m are the usual 
spherical harmonic coefficients in a geomagnetic model. 

For the satellite vector measurements, the field vector in earth fixed 
coordinates is assumed related to the orthogonal magnetometer axes, the 
bench coordinate system, fixed to the spacecraft boom by 

Bg = TgsTsmBm (27) 

Tqs is a function of the spacecraft attitude and position and relates the 
spacecraft coordinates to earth fixed coordinates. As presently 
configured, FIT can either assume this t ranf ormation is known or solve for 
its parameters. The assumed spacecraft coordinate system in FIT is such 
that the z axis is along the satellite velocity vector, the x axis is in 
the cross-track direction and the y axis completes the orthogonal system. 
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Note that this system differs from the APAFO system described elsewhere. 

The transformation Tqs is then calculated from the spacecraft inertial 
position and velocity, x$ and Jt s , and the Greenwich hour angle at the data 
times. Tsm is a function of Euler angles C defining the transformation 
from the orthogonal vector magnetometer axes to the spacecraft coordinate 
axes. The definition of the orientation parameters Ci, e 3 which 

describes the relationship between the magnetometer coordinate system and 
the spacecraft system are as follows: 

1) rotate by €2 about the z axis 

2) rotate by £3 about the new x axis 

3) rotate by about the new z axis 

The transformation matrix elements are then given by: 

(Tsm)11 = cos£icos£ 2 -cos£3sin£ 2 sin£i ; 

(T$m )12 = cos£isin£ 2 +cos£3COS£ 2 sin£i ; 

(Tsm) 13 35 sineisin£3 

(Tsm)21 = -sin£icos£ 2 -cos£3sin£ 2 cos£i • (28) 

(Tsm )22 * -sin£isin£ 2 +cos£3Cos£ 2 cos£i * 

(TSM ) 23 = coseisin£3; 

(Tsm) 3 1 = sine 3 sine 2 ; (Tsm) 32 * -sin£ 3 cos£ 2 ; (Tsm) 33 = cose 3 

To account for the possibility that the axes of the magnetometer are 
non-orthogonal , we relate the field vector projected on non-orthogonal unit 
vectors ui, u 2 , U3 to the orthogonal magnetometer axes. We assume that the 
ui vector is aligned along the x axis of the magnetometer, the u 2 vector 
lies in the x-y plane of the magnetometer system, and the U3 vector has 
arbitrary orientation . 

Let Bm represent the magnetic field vector in the orthogonal system and 
b N 0 = ( B l» B 2 1 B3) denote the projections of Bm along the ui, u 2 , U3 
vectors. Then 

Bl - Bmx 

B 2 = Bmx sina + BMy cosa ( 29 ) 

B3 = Bmx sin p COS7 + BMy sin fi sin 7 + Bm z cos P 

or Bno = TnmBM (30) 

where 

(TNM)ll = 1 ; (Tnm )12 = (Tnm )13 = ( t NM )23 = 0 ; 

(Tnm) 21 = Y i; (Tnm)22 = [I-Y1 2 ] 1 ' 2 ; 

(Tnm )31 = z i; ( t NM )32 = 2 2 (31) 

(Tnm )33 = [1 “ z l 2 ' Z2 2 ] 1 / 2 

and 


Yi = sin a-, Zi = sin/J COS7; Z2 = sin f) sin 7 


( 32 ) 
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The parameters Yi, Z \ , and Z 2 represent the non-orthogonal measurements to 
the orthogonal system, while the Euler angles €i, € 2 , e 3 relate the 


magnetometer system to the bench coordinates. Expressions for the partial 
derivatives used in FIT are found in the FIT documentation. 

The actual measurement of the magnetometer along an axis uj is a 
voltage, vj , and the field measurement, in nT, is 

B j = SjVj + (33) 

where the Sj and represent a slope parameter and bias, respectively. 

For convenience we choose to represent the voltage measurement in the form 

vj - fi jBj + flj (34) 

where fi j and jfj model the slope and bias parameters. Then 

▼ = TsiopeBNO + t* (35) 

The scalar magnetometer measurement is also represented with a bias 
parameter , 

v s = I Bg 1 + Tjs (36) 

The vector magnetometer measurement model is then: 

V = TsiopeTNM(TsM) T (TGS) T BG + fj . (37) 

Given a proper distribution of data, which takes one to two days to 


accumulate, and the assumptions mentioned above, the FIT program is then 
able to solve for the slope, non-orthogonality, bias offset and Euler angle 
parameters simultaneously with its solution for the spherical harmonic 
coefficients . 

In practice with Magsat, the slope and non-orthogonality parameters were 
equal to their nominal values, as expected, and were not solved for. 

However the Euler angles and bias parameters were an important part of the 
Magsat calibration. It is anticipated that this may also be the situation 
for APAFO. Phase one calibration should account for slope and non- 
orthogonality; indeed, finding fl = ^ = 7 = 0 in equations (29) would help 
confirm this (as well as optical bench and star sensor rigidity). 

If the spacecraft coordinate system used in T$m is fixed to the ESA/PPF 
spacecraft, then ( 61 , 62 ,^ 3 ) will be time-dependent due to boom flex and 
twist. Routine evaluation of these changing Euler angles may be a costly, 
and unecessary, source of attitude transfer error. If the boom indeed 
proves stable to within 0.5°, then it seems that dual star sensor 
determination of the orientation between the APAFO bench and the spacecraft 
will only be needed during calibration runs and, of course, independent 
boom stability checks. For routine evaluation of Tc,m and Tqs , the "S" 
coordinate system could then be fixed to the bench mounted APAFO star 
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sensors. Then TgM would represent the post-launch orientation of the bench 
system (fixed to the VFM optical cube) relative to the APAFO star sensors 
and will be steady to within the thermal stability of the bench systems. 
Given the position of a point on the spacecraft in geocentric coordinates 
and the values of Tqs from the APAFO star sensors, the position of the 
APAFO bench in geocentric coordinates can be determined to sufficient 
accuracy in the rigid boom approximation. 
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BOOH LENGTH /SPACECRAFT FIELD TRADE-OFFS 


Introduction 


The purpose of the boom is to place the magnetometer at sufficient 
distance from the spacecraft so that spacecraft fields are negligible. For 
APAFO this is not feasible. Rather, the boom length will be such that the 
spacecraft fields can be modeled with sufficient accuracy as to make them 
negligible. An error analysis study has been conducted to determine the 
trade-offs involved. 

This study is a preliminary error analysis for the APAFO magnetic field 
experiment. A vector magnetometer and a scalar magnetometer are located at 
the end of a boom. The spacecraft fields are due to currents in two torque 
rods and these currents are known at discrete time intervals. The 
positive x-axis is along the axis of the boom, which is perpendicular to 
the y-axis along the longitudinal axis of the spacecraft. The z-axis 
completes the right-handed reference system. The torquing rods are on the 
y-axis at (0.0, 7.2, 0.0) meters. One rod has its axis parallel to the x- 
axis . The other rod has its axis parallel to the z-axis. The length of the 
boom is expected to be 12 meters. See Figure 1. 

The total magnetic field is measured by a vector magnetometer and a 
scalar magnetometer at the end of the boom. When the boom is in its unbent 
nominal position, the vector magnetometer is located at (x, 0, 0) with axes 
parallel to the corresponding axes of the spacecraft reference frame. Star 
cameras are located on the spacecraft and at the end of the boom. At a 
given instant of time the ambient magnetic field is assumed to be constant 
in the neighborhood of the spacecraft; i.e., for an unbent boom the ambient 
field vector components are the same at the torque rods, the magnetometers 
and the star cameras. The maximum gradient of the Earth’s field is 
expected to be less than 28 nT/km, 24 nT/km in the horizontal plane. For a 
maximum distance of 20 m this gives a maximum field difference of 0.56 nT, 
0.48 nT in the horizontal plane. The ambient field difference across an 
optical bench of 2 m would be less than 0.06 nT. 

For this study, the spacecraft magnetic field is modeled by two magnetic 
dipoles, which have their origins at the centers of the torque rods. The 
field from these dipoles along the boom is shown in Figure 2. The dipole 
axes are coincident with rod axes and have corresponding polarities. 
Activation of the torquer rods will cause the boom to bend and oscillate 
relative to the spacecraft. Thus the position and orientation of the 
magnetometers and the outboard star cameras in the spacecraft coordinate 
system will change with time. In this study, the boom is assumed to have a 
parabolic shape in its fundamental vibrational mode and other modes are not 
considered . 

It is assumed that four independent magnetometer measurements are known 
at a fixed time: one scalar magnetometer measurement and three vector 
magnetometer component measurements. It is also assumed that apriori values 
are known for the parameters in the magnetic field models and the parabolic 
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boom model. Standard deviations are taken as known for the measurement and 
a priori errors. This information is used to compute a covariance matrix 
for the above parameters and the components of the ambient field using a 
weighted least-squares estimation scheme. 

The error sources considered in this study are 

1. magnetometer measurements (4), 

2 . torque rod magnetic dipole moments ( 2 ), 

3. orientation of the vector magnetometer relative to the spacecraft (3). 
The estimable parameters are 

1 . dipole moments for torquer rods ( 2 ), 

2 . orientation angles for the vector magnetometer relative to spacecraft 
(3), 

3. ambient magnetic field vector components (3). 

These parameters are not actually estimated in this study. Rather, the 
covariance matrix from the parameter estimation equations is computed using 
assigned values for the parameters and the standard deviations on the 
magnetometer measurements and the a priori parameter values. If the 
estimation process itself were performed, we would be solving in a weighted 
least-squares sense for eight parameters from nine data values (see error 
sources above). From the covariance matrix we obtain the standard estimates 
of error for the 8 parameters which would result if we actually solved for 
them. From this matrix we can also obtain the cross-correlations for the 
estimated errors in the parameters. 

We emphasize that the procedure described above is limited to data 
observed at a single fixed time. This is because models for the parameters 

as a continuous function of time are not assumed to be known. 

= x-component of the total magnetic field 

P 2 = y-component of the total magnetic field 

p 3 — z-component of the total magnetic field 

p 4 = magnitude of the total magnetic field. 

All the pi’s (i = 1 to 4) are in the vector magnetometer reference 
f rame . 

qi = magnetic dipole moment for torque rod parallel to spacecraft 
x-axis 

q 2 = magnetic dipole moment for torque rod parallel to spacecraft 
z-axis 

q 3 = a, the angle of rotation of the vector magnetometer about the 
x-axis of the spacecraft 

q4 = />. the angle of rotation of the vector magnetometer about the 
y-axis of the spacecraft 

q 5 = 7 , the angle of rotation of the vector magnetometer about the 
z-axis of the spacecraft 


32 


qg = x-component of the ambient magnetic field 

q 7 s y-component of the ambient magnetic field 

qg = z-component of the ambient magnetic field 

qg, q 7 and qg are in the vector magnetometer reference frame. 

In a solution for the eight parameters, allowing for a priori errors, a 
covariance matrix is obtained from which can be extracted the standard 
deviations of the estimated quantities. This covariance matrix is given by 


K = pTM-lp + N* 1 


where P T is the transpose of P. For the problem at hand P is a 4 x 8 matrix 
whose element in the i th row and column is 9pi/9qj evaluated for the 
nominal values of the parameters. The pj_ (i = 1 to 4) are the measured 
quantities. The qj (j = 1 to 8) are the model parameters: 2 torquer bar 
magnetic moments, 3 boom bending angles, and 3 ambient field vector 
components . 

M is a 4 x 4 diagonal matrix, the i th diagonal element being <7i 2 , the 

variance of the ith measurement (i = 1 to 4). N is an 8 x 8 diagonal 

matrix, the i^h diagonal element being s^Z, the apriori variance of the i th 
model parameter. 

Some of the questions which motivated this study are 

1. How accurately can the ambient field be determined? 

2. How does the ambient field accuracy vary with boom length? 

3. How does the ambient field accuracy vary with dipole model accuracy? 

4. What are the primary sources of error? 

The mathematical theory supporting this study is presented in the 
following sections. In the final sect ion, graphical results are given to 
support the answers given there to the foregoing questions. 


The Magnetic Dipole Model 

As a first approximation, the magnetic field produced by a torque rod may 
be represented by that of a magnetic dipole at the center of the rod and 
the positive dipole axis along the positive rod axis. The field of a dipole 
is obtained from the scalar magnetic potential 


V = m*u/4fu3 = m cos(5/4ju2, where (38) 

m = magnetic dipole moment vector, 
m = magnetic dipole moment (m = |m|), 

u = vector from the dipole origin to the point considered, 
u = magnitude of u (u = |u|), 

<5 = angle between the positive dipole axis and u„ 
m*u = scalar product of m and u. 
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The magnetic dipole field is 

B = grad v » where (39) 

fi 0 = permeability of free space* 

B can be broken into components 

B u = -f*o 8 V/ 9u = 2 B 0 cos<5, (AO) 

Bfl = -(/i Q /u) 9V/96 = B 0 sin$, where (Al) 

B 0 = /^o m /Airu 3 teslas* (A2) 


Since fl 0 = Ax x 10 ~ 7 Henrys /meter , we have 
B 0 = lO’^m/u 2 teslas. 

We convert this to nanoteslas by multiplying by 10 9 . Thus 
B 0 = 100 m/u^ nanoteslas. 

The units for m are ampere meters 2 (Am 2 ) if the magnetic dipole field is 
assumed to be produced by a single small closed loop of electric current. 
When the field is produced by multiple loops of current, it may be 
convenient to use ampere turns meters 2 (Atm 2 ) as units for m. 

x,y,z : coordinates of magnetometers in the spacecraft rectangular 
coordinate system 

V 1 > v 2 • v 3 : coordinates of a magnetic dipole center in the spacecraft 
rectangular coordinate system 

ui,U 2 ,U 3 : coordinates of the magnetometers relative to the dipole in the 
spacecraft rectangular coordinate system 


Ul = X 

- vi 

(A3) 

U 2 = y 

- v 2 

(AA) 

U 3 = z 

- v 3 

( A5 ) 

m^, m 2 , m 3 : components of magnetic dipole moment vector m in the 
rectangular coordinate system 

spacecraft 

m,0,X 

: spherical coordinates of m 


mi = m 

sin0 cosX 

( A6 ) 

m 2 - m 

sin0 sinX 

( A7 ) 

m 3 = m 

COS0 

( A8 ) 
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From (38) we get 


V = (miui + m 2 U 2 + 1113113 )/4tu3, (49) 

For the torque rod with an axis parallel to the x-axis in the spacecraft 

reference frame, we have from (46) through (49), 

V = (q 1 /4l')ui/u3. (50) 

The components of the magnetic field vector produced by this rod at the 
point (x, y t z) are then 

a X = -9V/9x = -100qi(l-3ui2/ u 2)/ u 3 nT, (51) 

a 2 - -9V/9y = 300qiuiU2/u^ nT, (52) 

33 = -9V/3z = 3 00 qiuiU 3 / u^ nT. (53) 

The computation of the covariance matrix requires values of the following 
partial derivatives . 

ci = 9ai/3qi = ai/qi (54) 

C 2 = 9a2/9qi = a 2 /qi (55) 

C 3 = 9a3/9qi = a 3 /qx (56) 

a ll = 9ai/9x = 300 qiui(3 - 5 ui 2 / u 2)/ u 5 (57) 

a 12 = 9ai/9y = 300 qiU2(l - 5ui^/u^)/u^ (58) 

ai3 = 9ai/9z = 300 qiU3(l - 5ui^/u2)/u5 (59) 

a2i = 9a2/9x =ai2 (50) 

a22 = 9a2/9y = 300 qiui(l - 5u2^/u^)/u5 (61) 

a 23 = 9a2/9z = - 1500qiuiU2U3 /u^ (62) 

a 31 = 9a3/9x = ai3 (63) 

332 = 9a3/9y = 323 (64) 

a33 = 9a3/9z = 300 qiui(l - 5u2^lu^-)lu^ (65) 

For the magnetic dipole with an axis parallel to the z-axi.s in the 

spacecraft reference frame, we have from (46) through (49), 

V = (q 2 /4lT)U3/u 3 . (66) 
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The components of the magnetic field vector produced by these dipoles at 
the point (x, y, z) are then 


bi = -9v/0x « 

= 300q2UiU3/u^ nT 

(67) 

b2 = -9V/0y > 

= 300q2U2U3/u 5 nT 

(68) 

b3 = -0V/0z < 

= -100q2(l-3u32/ u 2) / u 3 nT . 

(69) 

The computation of the covariance matrix will require values 
following partial derivatives. 

of the 

di = 9bi/9q2 

= t»i/q2 

(70) 

d£ = 9b2/9q2 

= b 2 / q2 

(71) 

d3 = 9b3/9q2 

= t»3 / q2 

(72) 

bn = 9bi/9x 

= 300 q2U3(l - 5 ui 2 /u 2 )/u5 

(73) 

bi2 = 9bi/9y 

= -1500q2UiU2U3 /u 7 

(74) 

b ^3 * 9bi/9z 

= 300 q2 u ld - 5u3 2 /u 2 )/u 5 

(75) 

b2i = 9b2/9x 

= bi2 

(76) 

b22 * 9b2/9y 

= 300 q2U3(l - 5u2^/u 2 )/u5 

(77) 

b23 = 9b2/9z 

= 300 q2U2(l - 5u3 2 /u 2 )/u5 

(78) 

b3i = 9b3/9x 

= t>13 

(79) 

b32 = 9b3/9y 

= b23 

(80) 

b33 = 9b3/9z 

= 300 q2U3(3 - 5u3 2 /u 2 )/u5 

(81) 

Let (hi, h£, h 3 ) = total spacecraft magnetic field in the 
reference frame and hij = Qhi/dxj. Then 

spacecraft 

hi = ai + bi 

(82) 

h2 = a£ + b£ 

(83) 

h3 = a3 + b3 

(84) 

hll = an + 

Ml 

(85) 

hl2 = a 12 + 

t>12 

(86) 

hi3 = ai3 + 

bi3 

(87) 
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( 88 ) 


h2l = 821 + b 21 


h 2 2 = 822 + b 22 

(89) 

h23 = 823 + b 23 

(90) 

h3i = 831 + b3i 

(91) 

h32 = 832 + b 32 

(92) 

b 33 = a 33 + b 33 • 

(93) 


A Model for the Boom Shape 

When torques are applied to the spacecraft, the boom will change shape 
deviating 0 from its static nominal position relative to the spacecraft The 
JooJoi the boom will remain fixed but the rest of the boom will oscillate 
about the nominal position. For this study, only the 

considered and the shape of the boom is assumed to be parabolic. The 
nominal position of the boom is along the positive x-axis of the spacecr 
reference frame with the foot of the boom at (0,0,0). The projection of the 
bo™ « an, given time into the xy-plane i. here deacrtbed by the equat.cn 


= ax 2 


(94) 


where a is a constant to be determined from the data. See Figure 3. 

Tf the tangent line at the end of the boom is projected into the xy- 
“ makesan angle 7 with the nominal boom axis and dy/dx = tan 7 . Let 

d = length of the boom and assume the projection of the magnetometers onto 
the *“Xs Is x = d, since the deviation of the magnetometers from the x- 
axis will be small. Values for 7 can be computed from star camera 
observations at the magnetometers and on the spacecraft. 


Differentiating (94) with respect to x gives 
dy/dx = 2ax. 

For small 7 , at the end of the boom x = d and dy/dx 

a = tan 7 /2d = 7 /2d. 

Substituting this in (94) gives 
y = 7 d/2 

at the end of the boom. 


(95) 

tan 7 = 7 . Thus 

(96) 

(97) 
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The projection of the 
equation 

z = bx 2 , 


boom into the xz-plane is here described by the 


(98) 


where b is a constant to be determined from the data. Using the same 
reasoning as before, with 7 replaced by /J, we find 

b = Pl2d, 

z = /7d/2 

at the end of the boom. 


(99) 

( 100 ) 


Measurements of p can be computed from star camera 
magnetometers and on the spacecraft* 


observations at the 


The boom may also twist about the x-axis through an angle a at the end 
of the boom. Measurements of a can be computed from star camera 
observations at the magnetometers and on the spacecraft. 


When the boom is bent from its nominal position, the values recorded by 
the magnetometers differ from those which would have been recorded in the 
nominal position. The y and z components of the magnetic field change due 
to the new position but the x component is, in effect, unchanged due to the 
smallness of the angles 7 and p. The perturbed position of the outboard 
magnetometers is (x,y,z) where 


x = d. 


( 101 ) 


y = 7 d /2 = q 5 d/ 2 , 


(102) 


z = P d /2 = q 4 d/ 2 . 


(103) 


The bending and twisting of the boom produces a reorientation of the 
vector magnetometer in the magnetic field. We assume, however, that the 
ambient field is estimated in the magnetometer coordinate system. This is 
equivalent to holding the magnetometers fixed in position and orientation 
and letting the spacecraft field sources move and twist. 


The vector magnetometer rotates about the x-axis through an angle a 
about the y-axis through an angle p, and about the z-axis through and angle 
7 : Since the angles are small, we can replace sina by a and cosa by 1 and 
similarly for p and 7 . 7 


We see now that the total magnetic field at the magnetometers is 

PI = 36 + hi - q 5 h£ + q 4 h 3 , 

P2 = q7 + h£ + qshi - q3h3 , 

P3 = 38 + h 3 - q^hi + q 3 h 2 , 


(104) 

(105) 

(106) 
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where a, /7, and 7 have been replaced by q3, q4 and qs. (hi,h2,h3) is the 
spacecraft magnetic field at (x,y,z) after translation [see (101), (102), 

(103)] but before rotation through angles (a,/f, 7) = (q3,q4*qs). ( qe * 
q8) is the ambient magnetic field vector and (pi, P2* P3) is the total 
magnetic field vector in the vector magnetometer reference frame. 

The computation of the covariance matrix requires values of the partial 
derivatives of the total magnetic field at the vector magnetometer [see 
(101) to (106)]. 

We differentiate pi, p2. and P3 in (104), (105), and (106) with respect 
to Q3» Q4» anc * Q5 an d then replace q3, q4, and q$ by their nominal values. 
For this study, the nominal values of q3, q4, and q5 are taken as 0, 
although in general, this will not be the case. For example 

9pi/9q4 = (9hi/0q4) - qs(9h2/9q4) + q4(9h3/9q4) + h 3 , 

0hi/9q4 = (0hi/9x) (0x/0q4) + (9hi/9y) ( 0y / 0q4 ) + (0hi/9z) (9z/9q4) = hi3d/2 
[See (101), (102), (103)], 

9pi /9q4 = hi3d/2 + h3 . 

The other partial derivatives are obtained in a similar way, giving 


P13 = 9pi/9q3 = 0, (107) 
P14 = 9pi / 0q4 = h]_3d/2 + h3, (108) 
P15 = 9pi/0qs = hi2d/2 - h2, (109) 
P23 = 9p2/9q3 = -h3, (110) 
P24 = 0P2/9q4 = h23d/2, (111) 
P25 = 0P2 /9q5 = h22d/2 + hi, (112) 
P33 = 0P3 /9q3 = ^2, (113) 
P34 = 9p3 / 9q4 = h33d/2 - hi, (114) 
P35 = 9p3 / 0q5 = h32d/2. (115) 


Obviously, pi6 = P27 = P38 = 1- 

The magnitude of the total magnetic field vector is 

P4 = (PI 2 + P2 2 + P3 2 ) 1 ^ 2 • (116) 
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If we let gi, g 2 , g 3 , and g 4 be the nominal values of pi, P 2 , P3. and P4* 
then 


(117) 

(118) 

(119) 

(120) 


gl = q 6 + hi, 
g 2 = q7 + h 2 , 
g 3 = q 8 + h 3 , 

g4 = (gl 2 + g2 2 + gl 2 ) 1 / 2 - 
From (116) we get the partial derivatives 

P 4 (0P4 /9qi) = Pl(0Pl/9qi> + P2(0P2/&qi) + P3(&P3/9qi)» d =1 to 8 * (121) 

With nominal values for the parameters, pi, (121) gives, for i = 3,4,5, 

P43 = (h 2 q 8 - h3q7)/g4, (122) 

P 44 = [(glhi 3 + g2^23 + g3 h 33 ) d / 2 +h 3q6 “ b iq8]/g4> (123) 

P45 = [(gl b 12 + g 2 b 22 + g3 b 32) d / 2 +b iq7 _ h2q6]/g4> (124) 

From (104), (82), (54), and (67), with nominal values for q^ and qs, we get 


Pll = ci = ai/qi . 

In a similar way, we get from (104), (105), and (106) 


P12 = d l = b i/q2. 

P21 = c 2 = a 2 /qi. 

P22 = d 2 = b2/q2. 

P31 = C 3 = a 3 /qi , 

P32 = d 3 = b 3 /q 2 . 

With nominal values for the parameters, 
P41 = (glPll + g2P21 + g3P3l)/84. 

P42 = ( glP12 + g 2 P 22 + g3P32 ) / g4 • 

From (104), (105), and (106) we get 


P16 = 1. 

P17 = 

P18 = 0. 

P26 = 

P27 = 1. 

P28 = 0. 

P36 = °* 

P37 = 0, 

P38 = 1- 


(126) 

(127) 

(128) 

(129) 

(130) 

pi, ( 121 ) gives, for i = 1 and 2 , 

(131) 

(132) 

(133) 

(134) 

(135) 
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With nominal values for the parameters, (121) gives, for i = 6, 7, and 8, 

P46 = gl/gA. P47 = g2 /gA • PA8 = g3/g4- (136) 

Results 


In the study, the spacecraft field was assumed to be due to two torquer 
rods with moments equal to 350 Am 2 each. The accuracy to which these 
moments are assumed known is Of. Of was allowed to assume the values 3.5, 
17.5 and 35 Am 2 or 1Z, 5Z and 10Z, respectively. The variation of Of, 
then, is a measure of how well the spacecraft field is known. Further, the 
angles between the torquer rods and the magnetometer platform were assumed 
to be known to an accuracy of o a , which for this study was assumed to be 
0.6° for a 12-meter boom. Both torquer bar moments and angles are assumed 
to be known a priori. For o& equal to 0.6°, the angular accuracy did not 
affect the error analysis. 

The measurement accuracies, denoted <7^ for the VFM and a sm for the SHM, 
were studied in the following combinations: (1) <7^ = 3.0 nT, o sm = 2.0 nT; 
(2) <?vm = 3.0 nT, <7 sm = 1.0 nT; (3) tfvm = 2 -° nT, <7sm = 1-0 nT. For each 
of these combinations, a set of calculations were performed and results 
showing the following were plotted: 

a) Uncertainty in the earth’s field vs. boom length for fixed 

Of (3.5, 17.5, 35 Am 2 ). Figures A, 5, 6; 1A, 15, 16; 2A, 25, 26. 

b) Uncertainty in the earth’s field vs. Of for fixed boom length 
(6, 12m). Figures 7,8; 17, 18; 27, 28. 

c) Uncertainty in spacecraft field vs. boom length for fixed 

Of (3.5, 17.5, 35 Am 2 ). Figures 9, 10, 11; 19, 20, 21; 29, 30, 31. 

d) Uncertainty in spacecraft field vs. Of for fixed boom length 
(6, 12m). Figures 12, 13; 22, 23; 32, 33. 

The x-component curve shapes are strongly influenced by the vanishing of 
the x-component of the input spacecraft field roughly 5 m out the boom. 

The y-component curve shapes are strongly influenced by the extrema in the 
y-component of the input spacecraft field roughly 4 m out the boom. 


Conclusions 


From these results it is possible to draw the following conclusions: 

!• Given minimal spacecraft field; i.e. long boom or very small Of, the 
accuracy of determination of Earth’s field is limited by the 
calibration accuracy of the VFM. The more accurately the VFM can be 
calibrated via the SHM, the more accurate the final measurements will 
be. (See, Figures A, 1A, and especially 2A.) 

2. There is a distinct trade-off between knowledge of the spacecraft 
field, Of, and boom length. If the spacecraft field were known 
perfectly no boom would be necessary. The actual trade-off depends, 
of course, on the actual strength of the spacecraft field. In this 
study, only the fields from the two torquer rods are included. 
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3. For the spacecraft field strength studied, if the field is known to 
5Z (17.4 Am 2 ), and if the VFM can be calibrated to 3.0-nT using SHM 
data good to but 2.0 nT, then a 12-m boom is long enough to meet the 
1 0-nT spacecraft field requirement (Figure 10). If the SHM data 
meets the 1.0-nT requirement, then a 12-m boom is long enough to meet 
some goals, but not all. It is long enough to meet the retirements. 
(Figure 25, also see Figure 15). 

4. For the spacecraft field strength studied, if the field is uown °. 

10Z (35 Am 2 ) a 12 _ m boom is not long enough to meet the goals is 

long enough to meet the requirement. (Figure 26, also see 16). 

5. If the spacecraft field is known to 1Z, then a 5-m boom length might 
be adequate. 


There will, of course, be other sources than the torquer rods. Some of 
those sources will probably not be anticipated prior to launch. The 
ability to accurately model the spacecraft field is the single most 
important factor in determining the necessary boom length. For known 
sources, this accuracy depends crucially on the adequacy of the model and on 
the accuracy and timeliness of the telemetry readout. The data rate 
required for the solid earth measurements is one per second. If telemetry 
readouts of the crucial spacecraft quantities are more widely spaced than 
this, or are not simultaneous with vector readouts, then some sort ot 
interpolation will be required which will reduce the accuracy. 

If sufficient attention is given pre-launch so that the model of the 
torquer rods, the solar panels, and any other comparable sources is 
certainly known and characterized by the available telemetry to better than 
1Z then the boom length could probably be shortened. The actual 
length would be subject to study but it is anticipated that it would not be 

shorter than 5 m. 
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MEETING THE ERROR BUDGET 


Two items from the error budget are pertinent, 
spacecraft field. 


SCALAR 


instrument error and 
VECTOR 



Req. 

Goal 

Req. 

Goal 

Instrument 

1.0 

0.7 

3.0 

1.5 

Spacecraft Field 

1.0 

0.5 

1.0 

0.5 

Rss 

1.41 

0.86 

3.16 

1.6 


Meeting the goals is feasible with the 12 m boom length, provided the 
vector instrument can be calibrated in-flight to 1.5 nT. Meeting the 
requirement should be possible provided the vector instrument can be 
calibrated in-flight to 3.0 nT. These both assume that the spacecraft 
field can be modeled to better than 52. As already noted, if a 12 model is 
possible, then the boom can probably be shortened. 


43 



APPENDIX A 


The MAGSAT tricoil system response was characterized by non-linear 
sensor cross-coupling, so Gjk* hence /Jjk> was a function of B and 12 
parameters hnm determined by ground calibration alone. The APAFO 4-coil 
system response is expected to have but weak linear cross-coupling because 
the field is to be nulled in each ring core plane rather than along but one 
axis. One might expect Gjk * ljk and /Jjk - a jk? however, the offset (and 
sense) coils 1 at x* 1 generate a stray sensor field B c (x*J,t) at sensor j t 
1; this field is distinguished from other ambient, spacecraft, and platform 
fields which comprise B(x*J,t) (due to proximity; the He scalar and the 
star cameras may generate a stray field at the fluxgate sensors and vice- 
versa which must meet mission requirements). So yj = Ekt a jk( B k + BC k)] P 
Ek[ajk B k]* But ® c depends upon the coil currents, so B c k = 

BC k(p n ( l^j ) , F(uii) , Fil) and the true field B + B c nulled by sensor j 
depends in part upon the coil currents, hence calibration, of the other 
sensors (Heisenberg strikes again) . To the extent that design ensures coil 
currents which are directly proportional to the true field being offset and 

| B c k I « I Bk I . B c k = El[bkl B l] --where Bi is assumed to be homogeneous 

across the fluxgates and bkl is determined by the geometric configuration 
of the fluxgates (notably the relative orientation and separation of the 
coils). Then yj = /Jjk B k = Ekt a jk( B k + Elbkl B l>] = Ek,l[ a jk B k + a jkt>kl B l). 
where /jjk is the response function operator. Reordering the summation and 
relabeling gives /Jjk = a jk + Eiajibik; thus the dependence of p on B can be 

ignored for an effectively linear response. Nevertheless, even if the 

response of sensor j were linear with no current in other coils, the 
presence of current in the other coils is expected to cause a response of 
sensor j . It follows that p should differ from an optical determination of 
a. Alternately, the effective sensor axis might be thought of as differing 
from the geometric axis of the coil; however, equations (5) for a need not 
apply to p and 12 independent elements of p, rather than 8 independent 
elements of a, need to be estimated. 

It may be possible to ensure that the dependence of /Jjk on B is 
negligible for the weak fields expected, but it seems wise to prepare for 
the contrary. Ground-based calibration must determine whether or not /Jjk 
depends upon B, and, if so, the functional form of such dependence-- 
including numerical values of h. Here it is assumed that /Jjk — /Jjk( B * a *h) 
with h given and /Jjk and B to be estimated iteratively. 
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